Topological phase transitions driven by gauge fields in an exactly solvable model 
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We demonstrate the existence of a new topologically ordered phase in Kitaev's honeycomb lattice 
model. This new phase appears due to the presence of a vortex lattice and it supports chiral Abelian 
anyons. We characterize the phase by its low-energy behavior that is described by a distinct number 
of Dirac fermions. We identify two physically distinct types of topological phase transitions and 
obtain analytically the critical behavior of the extended phase space. The Fermi surface evolution 
associated with the transitions is shown to be due to the Dirac fermions coupling to chiral gauge 
fields. Finally, we describe how the new phase can be understood in terms of interactions between 
the anyonic vortices. 
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I. INTRODUCTION 

Unlike other phases of matter, topologically ordered 
phases of many-body quantum systems can not be char- 
acterized by local symmetries. Instead, they are de- 
scribed by long-range properties, which are responsi- 
ble for the emergence of anyons, exotic quasiparticles 
with fractional statistics. Such phases occur in the cel- 
ebrated fractional Quantum Hall systems 1 or in a vari- 
ety of lattice models^. The latter offer microscopic con- 
trol and flexibility to study topologically ordered phases. 
Among them a pioneering role has been played by Ki- 
taev's honeycomb lattice model^, where the vortices can 
behave as non-Abelian anyons. The analytic tractability 
of the model has enabled the explicit demonstration of 
their various characteristics including fusion rules^, braid 
statistics^, topological degeneracy^, edge states^ and en- 
tanglement entropy-. Several experimental realizations 
of the model have also been proposed 9 -. 

Here we demonstrate the existence of a new topological 
phase in the honeycomb lattice model, which has been 
anticipated by KitaevJS. It appears for a fully packed 
lattice of vortices, when they behave as non-Abelian 
anyons. The ground state of this configuration is char- 
acterized by Chern number v — ±2, which implies the 
emergence of a new phase that supports chiral Abelian 
anyons. To understand the transitions between the dif- 
ferent topological phases, we study the evolution of the 
Fermi surface, which describes the long-range properties 
of the model. By considering the low-energy field the- 
ory of Dirac fermions, we show how two distinct types of 
topological phase transitions, i.e. different changes in the 
Fermi surface topology 11 , occur due to coupling to chiral 
gauge fields. We identify these to be due to dimerization 
and staggering of the model's couplings, both which can 
be realized in the proposed optical lattice experiments 9 . 
Finally, we derive analytically the phase boundaries for 
the new phase and illustrate how the transition between 
the non-Abelian and the chiral Abelian phase can be un- 
derstood in the context of anyon-anyon interaction driven 
phase transitions^. Due to an equivalence between the 
models^, our results apply directly also to p-w&ve super- 



FIG. 1: A brickwall representation of the bi-colorable honey- 
comb lattice, (a) The links are labeled x, y or z depending 
on their orientation, (b) Numbering of sites on plaquette p. 
(c) The elementary cell comprises of a black and a white site 
connected by a z-link. The cells are labeled by k = (m, n) in 
the orthonormal basis n x — (1,0) and n y — (0,1). (d) The 
full- vortex configuration is created by setting Jij — —1 on 
alternating z-links in x-direction (dashed links) and Jij = 1 
on all other links (solid links). The parallelogram shows the 
four sites belonging to the unit cell. 



conductors. 



II. THE HONEYCOMB LATTICE MODEL 

A. Mapping to free Majorana fermions 

Kitaev's model 3 , comprises of spin-1/2 particles resid- 
ing on the vertices of a honeycomb lattice. The spins 
interact according to the Hamiltonian 
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where Jij — (J a )ij are nearest neighbour couplings and 
a is either x, y or z when (ij) is an x, y or z-link, respec- 
tively (see Figure [Ha)). The second term is an effective 
magnetic field of magnitude AT, which explicitly breaks 
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time-reversal invariance. The sum in this term runs over 
all next to nearest neighbour triplets^. 

When the spin operators are represented as erf = ibfci, 
where Cj, b x , W and b\ are Majorana fermions, the Hamil- 

tonian takes the quadratic form H = | X)» j ^ij c i c j: 
where 

hij = 2J l jU ld + 2K)JjjkUjk, iii 3 = ib?bj. (2) 

The eigenstates | ty) of the original Hamiltonian ([J} are 
subject to the constraint 

A | *> = | *> , A = b^bVbla. (3) 

Since [H,Uij] = 0, the Hilbert space splits into sec- 
tors each labelled by u, a certain pattern of eigenvalues 
Uij = ±1. However, as A&y = — ?%A, should be 
understood as a classical gauge field with local gauge 
transformations A- The physical gauge- invariant sectors 
are labeled by the eigenvalues of the plaquette operators 
w p = U12U32W34M54M56M16, that satisfy [H, w p ] = (see 
Figure [ljb)). These can be identified with gauge invari- 
ant Wilson loop operators and their eigenvalues w p = — 1 
can be interpreted as having a vortex on plaquette p. We 
refer to the eigenvalue patterns u — {u^} and w = {w p } 
as gauge and vortex sectors, respectively. 

B. Gauge / coupling configuration equivalence 

To study the model as the vortex sector is varied, it 
is convenient to absorb the gauge u into the couplings. 
As can be seen from ([2]), the gauge choice — ±1 can 
be locally regarded as the sign of the physical coupling 
Jij. Therefore, the system can equivalently be viewed as 
being initially prepared in the vortex-free sector (w p — 1 
on all plaquettes) that contains the ground stated, with 
the other vortex sectors emerging by tuning the signs 
of the coupling configurations J = {Jij}. This view is 
also physically motivated as it relates the creation and 
transport of the vortices to the local manipulation of the 
couplings 5 . 

Strictly speaking one should also reverse the sign of 
the effecticve magnetic field couplings K locally. How- 
ever, when K approximates an external magnetic field, 
i.e. when K <C \J a \, considering only the signs of J a 
is a good approximation. Furthermore, we will see later 
that a non-zero K is required only to open a spectral 
gap. It does not affect the Fermi surface topology that 
characterizes the different phases. 

C. The full- vortex sector 

When the couplings satisfy J x m J y s» J z and K > 
and there are only a few vortices in the system, the vor- 
tices behave as strongly interacting non-Abelian Ising 
anyonsi The chiral Abelian phase emerges in this same 



coupling regime when one places a vortex on every pla- 
quette, i.e. considers the a full- vortex sector (w p = — 1 
on all plaquettes). It can be created, for instance, by 
using a gauge where Uij = ±1 alternates on all z-links in 
direction n-,,, while u^j = 1 on all other links. The lat- 
tice geometry and the corresponding unit cell containing 
four sites are illustrated in Figures [ljc)-(d). Normalizing 
I J a \ = 1, the gauge above is then equivalent to setting in- 
side the unit cell J x = 1, J y = 1 and (J z ,i, J z ,i) = (~ 1 5 l)j 
which corresponds to staggered couplings on z-links in di- 
rection n x . The system is translationally invariant with 
respect to {2xi Xl r\. v ) and a Fourier transformation gives 
the Hamiltonian H = j BZ d 2 p c^ p H p c pi where 

H P =( k h ? k %), (4) 

\ n bw rL bb / 

with 

abw ~ \ i i(_ e w- +e »p») J 

and 

h =K ( sm(p x -p y ) shx(p y ) ~icos(p x ) \ , g , 
bb \ sin(py) + icos(p x ) -sin(p x -p y ) J 

The integral is over the first Brillouin zone (— ^ < 
Px < f, — 7r < p y < 7r) and c p = 
(cb,p,Cb, P +irn lc ,Cw,p,c w ,p+irn K )- The indices (b,w) and 
(p> P + ^n^) refer to the two triangular sublattices of the 
honeycomb lattice and the two z-links inside the unit cell, 
respectively. 

This Hamiltonian can be readily diagonalizcd, with the 
lengthy expressions for the eigenvalues ±A,p and ±A,p 
given ins. In contrast to the vortex-free sector that has 
two Fermi points^, the low-energy behavior when K = 
is now dominated by the four Fermi points (A,Q = 0) 
located at = =p(f , f ) and Q| = ±(-f , ^). Their 
number and locations follow from the symmetries of Q, 

ri = <j y <E> 1 : r^pTl = h_ p , (7) 
r 2 = u z (g *y r 2 H p rl = h p+ ^ v , (8) 

which hold also for non-zero K. T\ acts only on the 
sublattice indices (b, w) and descibes a symmetry of the 
honeycomb lattice geometry. It is responsible for the 
doubling of Fermi points in the vortex-free sector. The 
new symmetry T 2 acts also on the indices (p, p + 7rn x ) 
that correspond to the two z-links inside the unit cell. 
Exchanging these links maps (J z , 1,^,2) = (— 1,1) — ► 
(1,-1), which preserves the full- vortex sector. Although 
T2 is a lattice symmetry, it can be equivalently viewed 
as an emergent global Z2 gauge symmetry as the corre- 
sponding gauges are inequivalent under the local gauge 
transformations A- This new symmetry is responsible 
for the further doubling of the Fermi points in the full- 
vortex sector. It changes the Fermi surface topology and 
thus implies a new phase. 
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FIG. 2: The normalized Hamiltonian @ defines a mapping 
from the Brillouin zone, which is topologically a torus T2, 
to a unit sphere S2. As K contributes to the low-energy 
theory only around the Fermi points, the map is reduced to 
the contributions from each. Depending on the sign of 
the term ±Ka z , the neighbourhoods of are mapped to 
either lower or upper hemisphere. As the ones with same 
orientations /j, end up to same hemispheres, the map winds 
around the unit sphere twice when viewed from the origin. 
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FIG. 3: A numerical study of the evolution of Fermi point 
locations (black dots) across the first Brillouin zone when (a) 
SJi = — > 00 (y = — 2 — > transition due to dimerization) 
and (b) <5J2 = 0— = — 2 — > — 1 transition due to 
removal of coupling staggering). The Fermi points move in 
the directions of the arrows and the annihilations occur when 
8 Ji = ^2 — 1 and 8 J| = \ , respectively. 



III. FERMI SURFACE ANALYSIS 

A. The low-energy theory of Dirac fermions 

To characterize this new phase and study its emergence 
from the phase with non-Abelian anyons, we consider 
its low-energy theory. The Hamiltonian Q is expanded 
around each Fermi point Q l ± by writing p = Q + k, with 
|k| <C 1. In general one obtains 

H Q =H Q + H%k x + H y Q k y + 0(k 2 ), (9) 

for some 4x4 matrices Hq. A projection onto the 2- 
dimensional low-energy subspace is given by 

Hq = PU Q H Q U Q P, (10) 

where Uq diagonalizes Hq and P projects onto the zero 

eigenvalue states of UqHqIIq. Keeping only the first or- 
der terms in k and projecting onto the low-energy states, 
we obtain 

i^^-k'T^^, (11) 

where o4 = {a x , ±(-1) V«), k l = { a %~^ v , ^"f" ) 

and cii = 2 — (— This Hamiltonian describes mas- 
sive Dirac fermions in two spatial dimensions. We in- 
terpret the mass being due to the scalar field K, which 
couples chirally, i.e. with a different sign, at the different 
Fermi points. 

B. Fermi surface topology and the Chern number 

The chiral coupling of K implies a chiral phase, which 
can be characterized by non-zero Chern numbe r 15 ' . We 
can relate the low-energy theory to the Chern number as 
follows. As the contribution of K vanishes away from 



the Fermi points, the four Hamiltonians (jll[) define an 
orientation preserving mapping from a torus (the first 
Brillouin zone) to a surface enclosing the origin (coordi- 
nates given in basis {cr a }). The degree of this map gives 
the Chern number— The orientation of Fermi point Q 
can be characterized by a winding number— 

VQ = ^-j> c Tr^iJpMi/p), (12) 

where Cq is a loop in the momentum space around Fermi 
point Q and T = a z (Si 1- Due to the chiral coupling of 
K, the neighborhoods of both (Q!_) with orienta- 
tions 1-j.Qi^ = +1 {[J-Qi = —1) are mapped to the lower 
(upper) hemisphere, with rest of the Brillouin zone being 
mapped to the equator. This is illustrated in Figure [2] 
For four Fermi points the map winds twice around the 
origin giving the Chern number v = — 2. We have verified 
this by calculating v also using the exact eigenstates^, 
as well as observed two non-degenerate edge states when 
a system of 10 3 lattice sites is placed on a cylinder—. The 
Chern number v = — 2 implies that the vortices in this 
phase are certain chiral Abelian anyons as cataloged for 
free fermion systems by KitaeA* 3 -. 

C. Emergent gauge fields and Fermi point 
transport 

The phase boundaries of this new chiral Abelian phase 
can be obtained by studying the effect of tuning the cou- 
plings J on the low-energy field theory. This allows one to 
obtain analytically the critical parameter values when the 
Fermi surface topology changes, i.e. when Fermi points 
are either created or annihilated. To this end we con- 
sider the system initially in the v = — 2 phase and set 
K = 0. The phase with non-chiral Abelian toric code 
anyons {y — 0) appears due to dimerization when, for 
instance, | J z \ >• \ J X \, \J y \ in Q 3 . This can be modelled 
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in the full-vortex Hamiltonian ((4]) as the perturbation 



8H X = -8J l0 - y ® a x 



(13) 



Projecting this into the low-energy subspace and combin- 
ing the linearized Hamiltonians for the paired Fermi 
points as -Hq; = diag(-Hq^ , HQi ), we obtain for small 

perturbations (8Ji <C 1) 



SHj = a 1 



(k 1 + 7 5 a;. 



(14) 



where now a 1 = (1 ® a x , (— 1) V z <g> a v ), 7 5 = a z ® 1 and 
= 8 J\ ^rxi (1; !)■ This Hamiltonian describes the 
Dirac fermions being coupled to a chiral gauge field 21 . 
As k* = (0, 0) no longer gives a vanishing Hamiltonian, 
the coupling to Aj shifts the paired Fermi points and 
towards each other, the direction being the same for 
both pairs. This agrees with 8H\ respecting both symme- 
tries ([7]) and (|8]), which also implies that all Fermi points 
have to vanish simultaneously. This is indeed the case as 
shown in Figure [3Ja,) , where we numerically demonstrate 
that dimerization in the large 8 J\ limit can cause local- 
ization of the fermions on the z-links and thus completely 
remove the Fermi points. 

We can similarly study the transition from the v = — 2 
phase to the non-Abelian Ising phase (y = —1). This oc- 
curs when one tunes the couplings between the staggered, 
{Jz,i,Jz,2) = (-1,1), and the uniform, (J z ,i, J z , 2 ) = 
(1,1), configurations. In (0} this is equivalent to the 
Hamiltonian perturbation 



5H 2 = 8J 2 (T y -t), 



(15) 



which respects the sublattice symmetry (|7|), but breaks 
©, the emergent symmetry responsible for the chiral 
Abelian phase. The low-energy theory is again a Dirac 
field coupled to a chiral gauge field (fl"4]) . but now A l 2 = 
^(l,(-l) l+1 V3), which shifts the pairs and 
independent of each other. This is confirmed by Figure 
&b), which shows that large 6J 2 distortions can cause the 
Qj_ Fermi points to annihilate while only transporting 
the other two. 

The study of the Fermi point transport holds also 
for K 7^ 0, which allows us to obtain the K depen- 
dent critical points 8Jf. As SJi is varied, the gapped 
Fermi points (the minima/maxima of the bands ±-Ei jP ) 
follow the same trajectories, although slower for larger 
K. The annihilations still occur at the critical momenta 
Q c = {(0,0), (f ,0), (f ,7r)}, where the gap always closes. 
This means that the eigenvalues of Hq c + SHj must 
satisfy ±Ei(K,SJj) = 0. The solutions are given by 

1 + 8J[ = V2 + X 2 and <SJ 2 C = which agree with 

our previous numerical studies^. In Figure 2] we outline 
the extended phase space as functions of 5J\ and 5J 2 
showing the three distinct topological phases. For gen- 
eral values of the couplings, the critical staggered config- 
uration and the tri-critical point occur when 
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FIG. 4: A section of the phase space as a function of SJi and 
8J2- The dashed lines are the phase boundaries and the circles 
are the locations of the K dependent critical points SJi = 
V2 + K 2 — 1 and 8 J% = . The phases characterized by 

the Chern number v = 0, v = —1 and v = —2 correspond to 
the dimerized, uniform and staggered coupling configurations, 
respectively. 



and 



Jt + K 



(17) 



respectively. These extend the previous results 2 ^ and 
show that a larger K has a stabilizing effect on the 
v = — 2 phase. 



IV. THE ROLE OF ANYON-ANYON 
INTERACTIONS 

In general, a transition from a non-Abelian Ising phase 
to a chiral Abelian phase has been predicted to arise due 
to anyon-anyon interaction s 12 ' 23 . As the Ising anyons in 
the honeycomb lattice model are strongly interacting 4 , 
we can explicitly establish this connection at an intuitive 
level. 

Based on numerical studies, we illustrate in Figures 
[SJa)-(d) how the spectrum evolves as the vortex den- 
sity inside a unit cell of fixed size is increased. Isolated 
vortices introduce modes with zero energy (Figure [SJa)), 
which acquire momentum dependance at close ranges due 
to interactions (Figure [5jb)). These zero modes describe 
the anyonic fusion degrees of freedom^, that are unique to 
non-Abelian anyons. When many vortices interact with 
each other simultaneously, these modes start forming a 
band structure (Figure [UJc)). Finally, as the density ap- 
proaches the limiting full- vortex sector, the characteristic 
four Fermi points emerge (Figure (5jd)). 

We observe that the Fermi surface is modified due to 
an emergence of a new low-energy band that is purely 
due to the interacting non-Abelian anyons. As demon- 
strated earlier, it is directly related to the Chern number 
that characterizes the ground state of the whole system. 
Therefore, our study reveals the microscopic mechanism 
for the anyon-anyon interaction driven phase transitions. 
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(a) (b) (c) (d) 

FIG. 5: A schematic illustration of the emergence of the 
full-vortex band structure due to interacting anyonic vortices 
(black squares) as the vortex density is increased, (a) A sepa- 
rated pair of vortices carries a zero mode (Ei :P — for all p). 
(b) Short-range interaction causes the zero mode to acquire 
momentum dependence, (c) The presence of many interact- 
ing vortices causes the zero modes to lose their character and 
to start forming a band structure, (d) The full-vortex band 
structure. 

V. CONCLUSIONS 

In summary, we have demonstrated the existence of a 
new topologically ordered phase in the honeycomb lat- 
tice model, which supports chiral Abelian anyons. This 
phase appears when the signs of the couplings alternate 
periodically, which is equivalent to creating a fully packed 
lattice of strongly interacting vortices. We characterized 
the phase by its distinct Fermi surface that arises due 
to an emergent lattice symmetry, and derived analyti- 



cally the phase boundaries by studying the Fermi point 
transport. We showed analytically that the transport, 
and thus the topological phase transitions, translate in 
the low-energy theory to coupling to chiral gauge fields. 
Their form is shown to be directly related to the lattice 
symmetries, which could be employed when engineering 
gauge fields in laboratory 24 , or studying the effects of 
lattice deformations in graphene 2 ^. 

The chiral Abelian phase can also be realized in the 
proposed optical lattice experiments, where the required 
staggered coupling configurations can be created by us- 
ing several parallel lattices with readily tunable spin 
interaction patterns 9 . Various experimental techniques 
can be considered for detecting the change in the Chern 
number—. Furthermore, as the v = — 1 phase with in- 
teracting non- Abelian anyons is equivalent to a p-wave 
superconductor—, a similar transition occurs also there 
for a sufficiently dense vortex lattice. Our results should 
apply as well to the chiral spin liquid variants of the 
model 2 ^. 
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